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ABSTRACT: Investigating gene expression of immune cells of whole blood or peripheral blood mononuclear cells (PBMC) under 
polyinosinic:polycytidylic acid (poly I:C) stimulation is valuable for understanding the immune response of organism to RNA viruses. 
Quantitative real-time PCR (qRT-PCR) is a standard method for quantification of gene expression studies. However, the reliability of 
qRT-PCR data critically depends on proper selection of reference genes. In the study, using two different analysis programs, geNorm and 
NormFinder, we systematically evaluated the gene expression stability of six candidate reference genes (GAPDH, ACTB, B2M, RPL4, 
TBP, and PPIA) in samples of whole blood and PBMC with or without poly I:C stimulation. Generally, the six candidate genes 
performed a similar trend of expression stability in the samples of whole blood and PBMC, but more stably expressed in whole blood 
than in PBMC. geNorm ranked B2M and PPIA as the best combination for gene expression normalization, while according to 
NormFinder, TBP was ranked as the most stable reference gene, followed by B2M and PPIA. Comprehensively considering the results 
from the two programs, we recommended using the geometric mean of the three genes, TBP, PPIA and B2M, to normalize the gene 
expression of whole blood and PBMC with poly I:C stimulation. Our study is the first detailed survey of the gene expression stability in 
whole blood and PBMC with or without poly I:C stimulation and should be helpful for investigating the molecular mechanism involved 
in porcine whole blood and PBMC in response to poly I:C stimulation. (Key Words: Pigs, Reference Genes, Whole Blood, Peripheral 
Blood Mononuclear Cells (PBMC), Polyinosinic:Polycytidylic Acid (Poly I:C)) 



INTRODUCTION 

Quantitative real-time PCR (qRT-PCR) can 
simultaneously measure gene expression in many different 
samples for a limited number of genes, and is especially 
suitable when only a small number of cells are available. It 
has become the standard method for quantification of gene 
expression studies. In case of qRT-PCR, the use of reference 
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genes as internal controls is the most common method for 
normalizing mRNA data (Huggett et al., 2005). The 
reference gene is a stably expressed gene that is 
experimentally verified in given species and tissues under 
given experimental conditions. However, an increasing 
number of reports have shown that the expression of some 
frequently used internal reference genes may vary 
significantly under different experimental conditions 
(Vandesompele et al., 2002; Nygard et al., 2007; Facci et al., 
2011). If the chosen reference gene has a large expression 
fluctuation among samples tested, the normalization will 
lead to erroneous gene expression profiles of the target gene 
of interest. Therefore, the selection of the most stable gene 
or set of genes as internal controls is a critical step to 
control the variability among samples for quantitative gene 
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expression studies with a sensitive qRT-PCR technique. 

Whole blood, especially peripheral blood mononuclear 
cells (PBMC) isolated from it, consists of many types of 
immune cells, such as lymphocytes, granulocyte and 
monocytes/macrophages, can produce different receptors 
and cytokines that lead to immediate innate immune 
responses (Gao et al., 2010; Uddin et al., 2012). As in vitro 
models, they have been widely used to dissect the 
pathogenesis and genetics behind infection and stimulation 
through gene expression studies (Gao et al., 2010; Uddin et 
al., 2012). Polyinosinic:polycytidylic acid (poly I:C) is an 
analogue of viral double-stranded RNA (dsRNA) that can 
serve as pathogen-associated molecular pattern and 
activates innate immunity through two major dsRNA 
sensors, toll-like receptor 3 (TLR3) and melanoma 
differentiation associated protein-5 (MDA5) (McCartney et 
al., 2009). Using the poly I:C as an immunologic stimulant 
is an important experimental measure to study the immune 
regulation of an organism to RNA viruses (Caskey et al., 
2011;Caoet al., 2013). 

Pigs are not only one of the most economically 
important livestock worldwide, but an important research 
animal that has become a useful model for various human 
diseases (Meurens et al., 2012). However, no studies have 
been performed to validate the stability of reference genes 
in samples of porcine whole blood and PBMC in response 
to poly I:C stimulation. Therefore, the aim of the current 
study was to identify a set of stably expressed reference 
genes in porcine whole blood and PBMC in response to 
poly I:C stimulation. The reference genes selected in the 
study would be helpful for investigating the molecular 
mechanisms involved in porcine blood and PBMC in 
response to poly I:C stimulation. 

MATERIAL AND METHODS 

Blood collection and PBMC isolation 

Three clinically healthy Landrace pigs of 35 days old 
were bled to obtain blood samples. For each piglet, 20 mL 
blood was collected via venipuncture into a vacutainer tube 
containing anticoagulant (EDTAK2). Out of the 20 mL 
whole blood per piglet, 15 mL was used for isolation of 
PBMC while the other 5 mL for whole blood culture. The 
whole study protocols for collection of the blood from 
experimental pigs were reviewed and approved by the 
Institutional Animal Care and Use Committee (IACUC) of 
Shandong Academy of Agricultural Sciences. 

PBMC were isolated using Ficoll-Hypaque PLUS (GE 
healthcare) following the manufacturer's instructions. In 
brief, the whole blood was firstly diluted by an equal 
volume of balanced salt solution (PBS). Then, 15 mL of 
diluted blood was carefully added on the top of 10 mL of 
Ficoll-Hypaque solution in a 50 mL conical tube and 



centrifuged at 460xg for 20 min at room temperature. After 
centrifugation, the middle whitish interface containing 
mononuclear cells was transferred to a new tube, and 
washed by adding 25 mL PBS followed by centrifugation at 
150xg for 10 min two times. 

Poly I:C stimulation 

PBMC, isolated from 15 mL blood of every pig, were 
diluted into 50 mL RPMI-1640 medium (Hyclone) 
supplemented with 10% fetal calf serum, streptomycin and 
penicillin, with the final cell concentration of ~2xl0 6 /mL. 
The cell suspension was divided into two parts, one of 10ml, 
used as the uncultured and un-stimulated control, the other 
40 mL, which was equally assigned into 16 wells (2.5 
mL/well) for the stimulation test. Then, these 16 wells were 
divided into four time-point groups, which were cultured 
for 3 h, 6 h, 12 h, and 24 h at 37°C with 5% C0 2 . 
Additionally, the four wells of every group were used for 
stimulation with different poly I:C concentration (0, 10, 20, 
40 |ig/mL, respectively) (Sigma-Aldrich). 

In order to assure a sufficient amount of cells for the 
next RNA extraction, samples of the three piglets with the 
same concentration of poly I:C and stimulation time were 
merged together for RNA extraction. Thus, we had 17 
samples of PBMC, including five samples without poly I:C 
stimulation (one uncultured control and 4 samples collected 
after 3 h, 6 h, 12 h, and 24 h of culture) and 12 samples 
with poly I:C stimulation (4 time points, and each point 
with 3 different poly I:C concentration tests). 

The other 5 mL fresh blood was diluted by adding 20 
mL RPMI-1640 medium supplemented with 10% fetal calf 
serum, streptomycin and penicillin. About 20 mL diluted 
whole blood was equally assigned into 16 wells (1.25 
mL/well) and the other 5 mL were used as the uncultured 
and un-stimulated control. According to the sample 
experimental treatments as PBMC, the diluted whole blood 
were stimulated and collected, and samples of the three pigs 
were merged. 

RNA extraction, cDNA synthesis and quantitative real- 
time PCR 

Total RNA was isolated from samples of whole blood 
and PBMC using the TRIzol reagent according to the 
manufacturer's protocol (Invitrogen, Burlington, Canada). 
The quantity and quality of RNA was determined by gel 
electrophoresis and NanoDrop 2000 spectrophotometer 
(Thermo Scientific). The PrimeScript RT reagent kit with 
gDNA Eraser (RR047A, Takara, Japan) was used to remove 
possible contaminating genomic DNA. The extracted RNA 
was firstly treated with gDNA Eraser for 2 min at 42°C, and 
then reverse-transcribed into cDNA for the mRNA 
reference gene analysis according to the instructions of the 
manufacturer. The final cDNA products were diluted 4-fold 
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prior to use in qRT-PCR. 

Six candidate reference genes, namely glyceraldehyde 
3-phosphate dehydrogenase (GAPDH), beta-actin (ACTB), 
beta-2-microglobulin (B2M), ribosomal protein L4 (RPL4), 
TATA-box binding protein (TBP) and cyclophilin A (PPIA), 
were selected from the most commonly used ones reported 
in previous studies (Martino et al., 2011; Cinar et al., 2012). 
The function, primers and other related information of these 
candidate reference genes are listed in Table 1 . 

All qRT-PCR were carried out using LightCycler 480 
SYBR Green I Master on Roche LightCycler 480 
instrument following the manufacturer's guidelines. Prior to 
performing the expression stability assay for all the primers 
of the six candidate genes, we carried out melting curve 
analyses to detect their specificity, and then calculated the 
PCR amplification efficiencies of every primer using a 
standard curve derived from a pooled cDNA mixture 
serially diluted 4-fold over five measuring points. In the 
expression stability assay, all the samples were amplified in 
duplicate. The PCR reaction consisted of 10 uL Blue- 
SYBR-Green mix, 1 uL forward and reverse primers (10 
pM/uL), 7 uL distilled water, and 1 uL of cDNA in a total 
volume of 20 |iL. The thermal cycling conditions were 3 
min at 95 °C, followed by 40 reaction cycles (10 s at 95 °C 
and 20 s at 60°C). The second derivative maximum 
algorithm included within the instrument software was used 
to determine cycle threshold (Ct) values for each reaction. 

Analyses of reference gene expression stability 

The data obtained from the qRT-PCR assays were 
converted into correct input files, according to the 
requirements of the software, and analyzed using geNorm 
version 3.5 (Vandesompele et al., 2002) and NormFinder 
(version 0.953) (Andersen et al., 2004). In brief, geNorm 
calculates the gene expression stability measure M for a 
reference gene as the average pairwise variation V for that 
gene with all other tested reference genes. Besides the M 
values, geNorm also calculates average expression stability 
values of remaining control genes during stepwise 
exclusion of the least stable control gene. Finally, to 



determine the optimal number of control genes for 
normalization, by the program, we calculated the pairwise 
variation V n /V n+1 between two sequential normalization 
factors NF n and NF n+ i containing an increasing number of 
reference genes. NormFinder uses the model-based strategy 
to estimate not only the overall expression variation of the 
candidate normalization genes, but also the variation 
between sample subgroups of the sample set, i.e. non- 
stimulation group and stimulation group in the study. The 
output of the program provides a stability value for each 
gene, which is a direct measure for the estimated expression 
variation enabling the user to evaluate the systematic error 
introduced when using the gene for normalization. 

Student's t-test analysis for unpaired samples was also 
performed to evaluate statistical differences between 
samples with and without poly I:C stimulation. 

RESULTS 

RNA quality and concentration of the PBMC and whole 
blood samples 

Total RNA of the samples of PBMC and whole blood 
was extracted using the TRIzol reagent, and their quality 
was measured by spectrophotometer. For the RNA of 
PBMC samples, the ratio A260/A280 was 1.94+0.25 
(average+SD), and the concentration 78.2 ng/uL+48.5 
(average+SD). For the RNA of whole blood samples, the 
ratio A260/A280 was 2.00±0.06 (average+SD), and the 
concentration 273.9 ng/uL+93.6 (average+SD). Our results 
indicated that all the RNA samples were good in quality and 
suitable for the next analysis. However, the concentration of 
RNA from PMBC was quite low, and larger volume 
samples should be used in the next cDNA synthesis. 

Verification of the primers for the candidate reference 
genes 

Using the mix pool RNA of the PBMC and whole blood 
samples, we generated amplification efficiency and melting 
curve analyses for the primers of the six candidate reference 
genes, and their results were shown in the last column of 



Table 1. Detail information of selected candidate reference genes for mRNA expression 



Gene 


Full name 


Function 


GenBank 
accession 
number 


Primer (5'-3') 


Reference 


Length 
(bp) 


qPCR 

efficiency 


GAPDH 


Glyceraldehyde-3-phosphate 
dehydrogenase 


Carbohydrate metabolism 


AF017079 


F: TCGGAGTGAACGGATTTG 
R: CCTGGAAGATGGTGATGG 


Martino et al., 
2011 


219 


97.90 


ACTB 


Beta-actin 


Involved in cell motility, structure 
and integrity 


DQ178130 


F: TCTGGCACCACACCTTCT 
R: GATCTGGGTCATCTTCTCAC 


Martino et al., 
2011 


114 


101.20 


B2M 


Beta-2-microglobulin 


Ctoskeletal protein involved 
in cell locomotion 


NM 213978 


F: TTCACACCGCTCCAGTAG 
R:CCAGATACATAGCAGTTCAGG 


Martino et al., 
2011 


166 


109.30 


TBP 


TATA binding protein 


Transcription initiation 
from RNA polymerase II promotor 


DQ178129 


F: GATGGACGTTCGGTTTAGG 
R: AGCAGCACAGTACGAGCAA 


Martino et al., 
2011 


124 


101.40 


RPL4 


Ribosomal protein L4 


Structural constituent of ribosome 


DQ845176 


F: AGGAGGCTGTTCTGCTTCTG 
R: TCCAGGGATGTTTCTGAAGG 


Cinar et al., 
2012 


185 


97.8 


PPIA 


Peptidylprolyl isomerase A 
(cyclophilin A) 


A member of the peptidyl-prolyl 
cis-trans isomerase (PPlase) family 


NM_214353 


F: CACAAACGGTTCCCAGTTTT 
R: TGTCCACAGTCAGCAATGGT 


Cinar et al., 
2012 


171 


95.3 
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Table 2. Expression levels of candidate reference genes for mRNA expression analysis with or without poly I: C stimulation 





H/XjJClllllCIlLal glOUp 


CA priH 


4 ("772 


DLiyi 


77? P 


RPT 4 

A / I.I 


PP1A 


PBMC 1 


Non-stimulation 


in i e+1 77 

JU.lO-Ll. / Z. 


90 7^+1 0f) a 


77 74.+1 27 










Poly I:C stimulation 


30.55±1.74 


22.1511.29" 


22.0711.18 


zy.jo±i.ij 


O ^ 1/14-1 £1 

Zj.l4±1.0z 


ZZ._o±l.UU 




All samples of PBMC 


30.4211.68 


21.6811.36 


22.1311.16 


29.3711.14 


25.3410.66 


22.3210.98 


Whole blood 


Non-stimulation 


31.5511.11 


22.0711.46 


24.2711.48 


31.3811.05 


28.6611.67 


24.5611.66 




Poly I: C stimulation 


30.9511.64 


22.4810.85 


23.5311.34 


31.4211.48 


28.8311.97 


24.1011.46 




All samples of 


31.1511.48 


22.3411.06 


23.7811.39 


31.4111.32 


28.7711.82 


24.2611.50 




whole blood 














All samples 




30.7911.60 


22.0111.25 


22.9511.52 


30.3911.60 


27.0512.45 


23.2911.59 


Expression levels of candidate reference genes were denoted by 


average+standard deviation (SD). 









Gene with different superscript indicates significant difference between gene expression levels in poly I:C stimulation group and non-stimulation group 
by t-test (p<0.05). 

' PBMC = Peripheral blood mononuclear cells. 



Table 1 and supplementary Figure SI (See e-version for 
supplement.), respectively. The amplification efficiencies 
for the six candidate reference genes ranged between 95.3% 
and 109.30%. The melting curves of the six pairs of primers 
were single, smooth and sharply defined with a narrow peak, 
suggesting all the six pairs of primers produce pure and 
homogeneous products. Thus, all the primers for the six 
candidate reference gene were specific and suitable for 
further qRT-PCR assays. 

Expression profiles of candidate reference genes 

The cycle threshold (Ct) values for the six (GAPDH, 
ACTB, B2M, TBP, RPL4, and PPIA) candidate reference 
genes are listed in Table 2. The Ct values for the six genes 
studied ranged from 20.75 (ACTB) to 31.55 (GAPDH). 
Compared between samples of whole blood and PBMC, the 
6 mRNA had similar expression trends. Out of the 6 mRNA 
candidate reference genes, ACTB (Ct 22.01), B2M (Ct 
22.95) and PPIA (Ct 23.29) had high expression levels, 
followed by RPL4 (Ct 25.34), TBP (Ct 29.37), and GAPDH 
(Ct 30.42) with the lowest level. When expression values 
were compared between samples with or without Poly I:C 
stimulation, only ACTB was expressed differentially (p< 
0.05) in PBMC samples, while no significant difference was 
observed for the other genes in the whole blood as well as 
in PBMC samples. 

Gene expression stability assessment using geNorm 

Gene expression stability measure M was calculated for 
the samples of whole blood and PBMC with or without 
poly LC stimulation. Genes with the lowest M value have 
the most stable expression, while the highest M value has 
the least stable expression. The expression M values of the 
six reference genes are plotted in Figure 1. It can be seen 
that, except more stably expressed in samples of whole 
blood than those of PBMC, the six candidate reference 
genes perform similar trend of expression stability in the 
sample of whole blood and PBMC. That is, B2M is the most 



stable reference gene with the lowest M value (0.631 and 
0.777 in whole blood and PBMC, respectively), while 
ACTB the least stable one with the highest M value (0.921 
and 1.426 in whole blood and PBMC, respectively). 

Assuming that the reference genes are not co-regulated, 
stepwise exclusion of the gene with the highest M value 
results in a combination of two constitutively expressed 
reference genes that have the most stable expression in the 
tested samples (Vandesompele et al., 2002). Figure 2a and 
2b shows the average expression stability values of 
remaining control genes for whole blood and PBMC. 
Combining using the two reference genes for normalization, 
B2M and PPIA, the average expression stability values can 
obtain rather small M values for the samples of whole blood 
(0.299) and PBMC (0.490). 

In addition to a ranking regarding the expression 
stability of each gene, geNorm also calculates a 
normalization factor for each sample based on the most 
stable reference genes, which are illustrated in Figure 3. A 
large variation means that the added gene has a significant 
effect on the normalization and should preferably be 
included for calculation of a reliable normalization factor. 




0.2 ■ 

0 -I 1 1 1 1 1 1 

ACTB RPL4 GAPDH TBP PPIA B2M 

Figure 1. Gene expression stability (M values) of the six 
candidate reference genes calculated by geNorm. 
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Average expression stability values of remaining control genes 
for whole blood 




Average expression stability values of remaining control genes 
for PBMC 



1 - 




0 4 1 1 ' ' 1 

ACTB RPL4 GAPDH TBP B2M 

PPIA 

<::::: Least stable genes Most stable genes ::::> 



(b) 

Figure 2. Average expression stability values (M) of remaining control genes during stepwise exclusion of the least stable control gene, 
(a) whole blood; (b) peripheral blood mononuclear cells (PBMC). 



We observed the lowest variation V„/V n+1 between inclusion 
of the fourth and fifth most stable reference gene (0.136 and 
0.149 for the whole blood and PBMC). 

Gene expression stability assessment using NormFinder 

The gene expression stability was also analyzed by 
another analysis program, NormFinder, where stable gene 
expression is indicated by lower average expression 
stability values. According to NormFinder results as shown 
in Figure 4, the six candidate genes were more stably 
expressed in whole blood than PBMC, and also had a 
similar trend of expression stability in the sample of whole 
blood and PBMC. Especially, TBP was ranked as the most 
stable reference gene in both whole blood and PBMC, 
while B2M and PPIA were the second and third most stable 
ones. NormFinder found the stability value for best 



combination of two genes, TBP and PPIA, were 0.636 and 
1.623 for whole blood and PBMC respectively. In 
accordance with results valuated by geNorm, the frequently 
used reference genes ACTB and GAPDH exhibited 
substantial variation, and were among the three least stable 
genes. 

DISSUSSION 

Whole blood and PBMC, consisting of many types of 
immune cells are important in vitro models to investigate 
the immune response to stimulation by environment or 
pathogens. As an analogue of viral double-stranded RNA, 
poly I:C is widely used as immunologic stimulant to study 
the immune regulation of an organism to the RNA viruses. 

Previous studies have indicated that different reference 
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Determination of the optimal number of control genes for normalization 
for whole blood 
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Determination of the optimal number of control genes for normalization 
for PBMC 



s 0.200 

o 
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S 0.150 

<u 
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(b) 

Figure 3. Pairwise variation (V„/V ll+1 ) analysis between the normalization factors NF„ and NF n+i . (a) whole blood; (b) peripheral blood 
mononuclear cells (PBMC). 



genes were required for gene expression studies in PBMC 
at different times or with different immunologic stimulants. 
For instance, RPL4 and TBP ranked as the two most stably 
expressed genes in PBMC of pigs at day 1, day 7, day 14 
after heat stress, while B2M was the most stable at day 21 
after heat stress (Xiang-Hong et al., 2011); PPIA, B2M, and 
RPL4 were the most stable reference genes in 
Lipopolysaccharide (LPS) stimulated PBMC, while 
Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase 
activation protein, zeta polypeptide (YWHAZ), RPL4 and 
PPIA were the most stably expressed reference genes in the 
case of lipoteichoicacid (LTA) stimulated PBMC (Cinar et 
al., 2012; Cinar et al., 2013). So far, no studies were 
performed to validate the stability of reference genes in 
cultured porcine whole blood and PBMC in response to 
poly I:C stimulation. 

There are a number of programs that allow assessment 



of multiple reference genes, such as geNorm 
(Vandesompele et al., 2002), NormFinder (Andersen et al., 
2004) and Bestkeeper (Pfaffl et al., 2004). But to date, there 
is no consensus on which method should be used to 
examine reference gene expression stability. In particular, 
Radoni et al. (2005) suggested that ranking lists of the 
genes tested was tool-dependent. A comparison of different 
algorithms of reference gene selection allows a better 
evaluation of the most reliable controls and reduces the risk 
of artificial selection of co-regulated transcripts. Therefore, 
in the present study, we used two statistical approaches, 
geNorm and NormFinder, to evaluate the gene expression 
stability of the six reference genes, and also compare the 
results obtained using these two programs. 

Based on the geNorm programs, the most stable 
reference genes in both whole blood and PBMC were PPIA 
and B2M. Combination of the two genes can obtain rather 
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ACTB RPL4 GAPDH PPIA B2M TBP 



Figure 4. Gene expression stability of the six candidate reference 
genes calculated by NormFinder. 

small M values for the samples of whole blood (0.299) and 
PBMC (0.490). As experimental data from the geNorm 
developers have shown that stably expressed genes 
typically exhibit mean M values lower than 0.5 in relatively 
homogeneous sample panels (Hellemans et al., 2007). Thus, 
the combining of the two genes, PPIA and B2M, is 
appropriate for normalization of gene expression of whole 
blood and PBMC. On the other hand, based on the 
NormFinder, TBP was ranked as the most stable reference 
gene in both whole blood and PBMC, while genes B2M and 
PPIA were the second and third ones. This discrepancy is 
probably due to the different algorithms used in the two 
programs. Though the most stable genes were different, the 
first three most stable genes were the same as evaluated by 
the two programs. 

It has been suggested that using the geometric mean of 
two or more reference genes for qRT-PCR studies might 
generate more reliable results (Vandesompele et al., 2002). 
In order to determine the optimal number of required 
reference genes, we also performed pairwise variation 
(V n /V n+ i) analysis by geNorm, as shown in Figure 3. Based 
on the cut-off value of 0.15 proposed by geNorm program, 
four reference genes are necessary to obtain the lowest 
control genes in the tested samples of both whole blood and 
PBMC. In practice, it is not always possible to measure 
multiple reference genes due to limited sample availability 
and cost. Additionally, the cut-off value of 0.15 is an 
empirical value and should not be taken as a too strict cut 
off value, as it is already suggested by the geNorm manual 
itself. Therefore, we recommend using the geometric mean 
of the first three stable genes, TBP, PPIA, and B2M, to 
normalize the gene expression of whole blood and PBMC 
with poly I:C stimulation. 

Up to now, most of the published qRT-PCR results were 
normalized to either GAPDH or ACTB without 



experimentally confirming their invariant expression under 
the experimental conditions described. In our present 
analyses, these two genes were the least and the third least 
stable genes among the six candidates. Consistent with our 
result, some previous studies also reported ACTB and 
GAPDH were unsuitable as reference genes in whole blood 
(Brym et al., 2013) and PBMC (Spalenza et al., 2011; 
Xiang-Hong et al., 2011) and other specific samples (Uddin 
et al., 2011; Zhang et al., 2012). ACTB is one of the major 
components of cytoplasmic microfilaments in eukaryotic 
cells and plays an important role in diverse cellular 
functions, while GAPDH not only acts as a component of 
the glycolytic pathway but also takes part in other processes 
as well. Both ACTB and GAPDH have sophisticated 
functions in organisms, and their expression profile might 
fluctuate according to the corresponding experimental 
conditions. Our analysis results confirmed, once more, 
blindly choosing reference genes is no longer acceptable, 
and authors must accurately validate candidate reference 
genes before use in gene expression qRT-PCR studies. 
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